A Constrained-Path Quantum Monte-Carlo Approach 
for the Nuclear Shell Model 
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A new Quantum Monte-Carlo (QMC) approach is proposed to investigate low-lying states of nuclei 
within the shell model. The formahsm rehes on a variational symmetry-restored wave-function to 
guide the underlying Brownian motion. Sign/phase problems that usually plague QMC fermionic 
simulations are controlled by constraining stochastic paths through a fixed-node like approximation. 
Exploratory results in the sd and pf valence spaces with realistic effective interactions are presented. 
They prove the ability of the scheme to yield nearly exact yrast spectroscopies for both even- and 
odd-mass nuclei. 
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The shell model with configuration interaction [1] is 
one of the most powerful frameworks for nuclear struc- 
ture calculations. In this picture, the nucleons beyond an 
inert magic core are confined in an active shell and in- 
teract through an effective two-body residual potential. 
The wave-functions of the nucleus are then determined 
by diagonalization of the Hamiltonian in the set of all 
the possible configurations and used to calculate physi- 
cal observables. Unfortunately, the exponential scaling 
of the many-body space with the number of valence nu- 
cleons or the size of the single-particle basis puts strong 
restrictions on the applicability of the shell model. 

Quantum Montc-Carlo (QMC) methods potentially of- 
fer attractive alternatives to such limitations. The many- 
body problem is indeed reduced to a set of stochas- 
tic one-body problems, numerically tractable, describing 
independent particles that randomly walk in fluctuat- 
ing external fields. In this way, exact correlated wave- 
functions are reconstructed by averaging independent- 
particle states called walkers. To date, the Shell Model 
Monte-Carlo method (SMMC) is the main application of 
a QMC approach for the nuclear shell model [2]. It re- 
lies on the standard sampling by auxiliary-fields of the 
path-integrals resulting from the Hubbard-Stratonovich 
transformation of the imaginary-time propagator. With 
schematic residual interactions, the SMMC method ac- 
curately reproduces the properties of nuclei at zero and 
finite temperature [3]. However, with realistic effective 
interactions, the approach is plagued by the so-called 
fermion sign/phase problem which signature is an expo- 
nentially vanishing signal-to-noise ratio as the tempera- 
ture is decreased. To overcome this pathology, a family of 
modified Hamiltonians is usually sampled where interac- 
tions leading to the problem are artificially reduced. The 
SMMC results have then to be extrapolated to recover 
exact observables [4]. Moreover, the SMMC method may 
not be used to obtain a detailed spectroscopy of nuclei. In 
contrast, low- lying states can be reconstructed as long as 
the stochastic process is used only to generate a subspace 



in which the Hamiltonian is diagonalized. This so-called 
Monte-Carlo Shell Model technique [5] has been success- 
fully applied to very large configuration spaces and pro- 
vides variational estimates of nuclei properties at low en- 
ergy [6]. 

In this Letter, we present a new QMC scheme to inves- 
tigate yrast states in the framework of the shell model. 
The sign/phase problem is managed via a constraint on 
the Brownian motion involving a trial wave- function with 
restored symmetries in the spirit of fixed-node ab initio 
calculations [7] where the random walk takes place in the 
position space. 

Let us first consider a generic two-body Hamiltonian 
H cast in a quadratic form of one-body operators T and 
{Os}: 
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z] (Cj) is the creation (annihilation) operator of 
a fermion in a single-particle state \i) of an orthonormal 
basis. This decomposition can be achieved by performing 
a Pandya transformation, as detailed in Ref. [2]. 

All QMC approaches use a stochastic representa- 
tion of the imaginary-time propagator exp(— rJ?) to 
project an initial wave-function |$o)i supposed here to 
be a Slater determinant, onto the ground state |^'g): 
limT-_>oo exp(— r_ff)|$o) I^'g)- In order to improve the 
efficiency, Zhang and Krakauer [8] have proposed to bor- 
row the idea of importance sampling by generating walk- 
ers |<I>) according to their overlap with a trial state |^'t) 
not orthogonal to the ground state. By directly including 
the complex importance function (\E't|$) in the Brown- 
ian motion of walkers, one is led to reformulate the exact 
state at a given imaginary time r according to 



exp 



E 



1$. 



(2) 



where E[-] denotes the average of a stochastic process. 
The weight H^. is a Boltzmann-type exponential factor. 
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and the evolution of the single-particle states {|<^n)} of 
the Slater determinant |<l>) are driven by the following 
Ito stochastic diflFerential equations: 



= (*T|*o)exp 



d|</-„) = [-dr(T-J2uj,[0',+{6,)^^^^Os]) 



(3) 



\<Pn). 



Here, we introduced the so-called local estimate 
(i)^^^ = {^t\A\^)/{^t\^) of any observable A. The 
{VKjj} refer to independent Wiener processes with vanish- 
ing ensemble averages and the following multiplication 
table for their infinitesimal increments [9] 



Vs,s' : E[dWs 
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The stochastic dynamics defined by (3) differs from the 
underlying one of the standard SMMC method by the lo- 
cal estimate of the Og operators in the drift term. There- 
fore, through this additional term, guides the ran- 
dom walk towards a region of the over-complete manifold 
of Slater determinants where the contribution to |^g) is 
expected to be large. 

From a general point of view, any QMC sampling 
is plagued by sign/phase problems when the phase of 
the overlap varies during the Brownian motion. 

In that case, a population can yield a mean overlap 
E[(^'g|$)] = 0, the contributions of these realizations 
cancelling each others. Hence, such trajectories are use- 
less for reconstructing the ground-state wave-function, 
and they severely degrade the signal-to-noise ratio be- 
cause they only increase the statistical error. Within the 
QMC scheme (2), the overlap is given by 

(*g|$> cx (*T|exp(-Ti?)|$) = E[n,], (5) 

Ilr being defined by (3) but with Hq = (5't|$)- Control- 
ling the phase problem requires to use an approximation 
that will constrain the random walk by a change E ^ 11 
of the weight of the realizations. First, by taking the 
real part of the local energy in its evolution (3), we en- 
sure that the overlaps with the ground state and 
with the trial state have the same phase. Now, we 
have to prevent the random walk from populating sym- 
metrically the complex (4'T|$)-plane to avoid the mean 
overlap merging with the origin. To this end, we apply 
the phaseless approximation [8], where the constrained 
weight n is supposed to evolve according to the dephas- 
ing dd = a.Tg{^T\^T+dT) / {^t\^t) of the overlap with 
the trial state as 

IlT+dT = Ht exp^— dTRe(J?)^^^ ^ max{0;cos(d(?)}. 

(6) 



Therefore, the more the phase of changes during 

dr, the more the weight of the associated walker is de- 
creased, and when |d^| > 7r/2, the walker is discarded. 
In this way, the centroid of the population is maintained 
into the half- plane Re(^'T|4') > even if some realiza- 
tions having reduced weights are generated in the other 
half-plane. This constraint thus offers a compromise be- 
tween the need to control the phase problem and the 
conservation of the initial form for t^he probability dis- 
tribution. Finally, replacing H by H in (2) allows one 
to reconstruct jvithout phase problem an approximate 
ground state \^g)- We emphasize that approximation 
(6) also removes the risk of an infinite variance on the er- 
ror I'^g) ~ I^)) which would result from an accumulation 
of realizations close to {^t\^) = 0. 

The phaseless QMC method described above has been 
applied to many-body problems in quantum chemistry 
[8], with a simple mean-field wave- function for the trial 
state I^'t). For the nuclear shell model, we propose here 
an improved scheme to reconstruct the ground state in 
each angular-momentum channel J, M. Starting from a 
Slater determinant |$o)i ^ satisfying test wave-function 
to initiate, guide, and constrain the dynamics can 
be obtained by mixing all the possible values of the spin 
projection K in the intrinsic frame: 



.7 

K=-J 



ciPi 



(7) 



Here, the operators P'mx '"^^^ weighted averages of rota- 
tions Uq. parametrized with Euler's angles Q. as [10] 



pJ _ 



2J-H 
87r2 



(8) 



where Dj^j^ denotes Wigner's D-function. Below, we re- 
strict the Slater determinant |<I>o) to be factorizcd into 
a product of independent-neutrons and -protons wave- 
functions so that the ansatz (7) has also a good isospin 
projection. Moreover, when the valence space contains 
only one major shell of the harmonic oscillator, the par- 
ity is already a good quantum number and no further 
restoration is needed. Rotations Ua simply transform 
|$o) into another Slater determinant. As a superposition 
of symmetry-related Slater determinants, l^fj.*^) is no 
longer an independent-particle state and can absorb cor- 
relations between particles beyond the mean-field level. 
The energy Eq of an yrast state I^'g*^) for the angular 
momentum J can now be determined after a large enough 
imaginary-time as 



Eq - 



JM\ 



(^j,|ge-"^|$o) 



(9) 



provided the approximate state j^Pj,*^) and the exact 
wave- function have a non-zero overlap. In Eq. (9), 
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we have noted 1*^) = Y.K.K' <^'k'^k'Pkk'\^) ^he 
new state that emerges from rotational invariance and 
from the well-known properties {Pmk)^ 



pj p.T' 



Prm and 



= ^jj'Skk'P'mm' ■ This wave- function ap- 
pears as a natural candidate to apply the QMC formalism 
(2)-(3) with the biased weight (6) to constrain stochastic 
paths. An approximate yrast state for any spin J 

is thus obtained with energy Eq deduced from (2) and 
(9) 



En 



E 



(10) 



Note that a similar form also holds for any scalar obscrv- 
ables commuting with the Hamiltonian. In other cases, 
the well-known mixed estimate ^ gives an ap- 

proximate ground-state expectation value of an observ- 
able A: 



, V. (mix)_ 



Re(^^^|A|^^^) 
Re(*^^|*^^) 



(11) 



It can be corrected by the extrapolated estimator 



magnitude better in the difference - \^t^). 

In zero-temperature QMC methods, the imaginary 
time needed to converge to the ground-state wave- 
function depends on the quality of the approximate state 
used to initiate the propagation. In the phaseless QMC 
scheme, this state also guides and constrains the 

motion of walkers. All these considerations naturally 
lead to the choice of the variational solution obtained 
by energy minimization in the subspace of Slater de- 
terminants after quantum number projections. Such a 
symmetry-entangled mean-field (SEMF) method was re- 
cently applied to the two-dimensional Hubbard model 
in condensed matter physics [11] and in quantum chem- 
istry [12]. In nuclear physics, the SEMF method is sim- 
ilar to the so-called VAMPIR approach [13], but with a 
Slater determinant instead of a quasi-particle vacuum, 
and thoTigh the energy minimization is carried out in a 
different way. First, we note that the energy E:^ in the 
trial state (7) reads 



that is one order of 



(12) 



so that variation with respect to the amplitudes C^* 
yields a generalized eigenvalue equation: 



K' 



K' 



(13) 

On the other hand, owing to Wick's theorem [14], the 
expectation values of HP^j^, and P^k' |$o) are only 



expressed with the contractions {po)ij — (c^Cj)^ , po be- 
ing the one-body density matrix associated with |^'o)- 
Therefore, the energy (12) is a functional of po, i-e. 
Ej, = Ej,[pq]. The stationarity condition under the con- 
straint of orthonormal single-particle states {|^o,n)} im- 
mediately leads to the following Hartree-Fock like equa- 
tion 



/i-^[po],Po =0 



(14) 



The effective single-particle SEMF Hamiltonian is de- 
fined as the gradient of the projected energy according 
to the density, h'fj = dE;^/d{po)ji- This derivative can 
be further calculated by using Wick's theorem extended 
for matrix elements [14] through the introduction of one- 
body density matrices TZn between |<I>o) and C/n|$o)- The 
SEMF Hamiltonian is finally given by: 



h^ipo] = JdnxJ^A^^ [{Un - 1) {£[nn] - E^) 



EKK'C'K*C'K'Di*K'mdet{An) 



(15) 



Here, Aq_ and _Bo depends on po and on the matrix rep- 
resentation Uq, of rotations in the one-particle space ac- 
cording to Aq = l + {Un-l)po and En = l-l-po(C^n-l); 
5[7?.o] and h\R.n\ stands for the usual Hartree-Fock en- 
ergy functional and Hamiltonian, but in which the den- 
sity is replaced by TZn- This transition density matrix 
can be expressed as TZq = poC^o^fj"'^- More details about 
the derivation of (15) can be found in Refs. [12, 15]. 

Let us now present phaseless QMC results for the nu- 
clear shell model. We address here systems for which 
exact results are available for benchmark. Fig. 1 displays 
the SEMF, QMC, and exact yrast spectra of three kinds 
of s(i-shell nuclei with the USD interaction [16]: K N = Z 
odd-odd nucleus (^®A1), an odd- mass nucleus (■^''Na), and 
an even-even nucleus (^^Mg). Also shown is the spec- 
trum of ■''^Ni in the pf model space with the GXPFIA 
interaction [17]. Exact energies are extracted from [18]. 
First, we observe in Fig. 1 that the SEMF method of- 
fers a good approximation in all considered cases: Yrast 
energies differ from the exact values by less than about 
1 MeV. Second, relative errors on the QMC energies at 
any angular momentum do not exceed 0.3 % with statisti- 
cal error bars about 40-50 keV. In particular, the results 
for ^''Na give good evidence that the phase problem is 
well controlled in our calculations. Indeed, the odd-mass 
nuclei are more pathological with respect to QMC simu- 
lations, because of the time-reversal symmetry breaking 
of the stochastic dynamics that leads to a phase problem 
irrespective of the interaction. 

To further test the nuclear phaseless QMC scheme, we 
now examine the neutron (proton) occupation numbers 
iV"(p) of the various shells (nlj), the electric quadrupole 
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FIG. 1. (Color online) Yrast spectra of three sd-shell nuclei 
and one p/-shell nucleus as obtained with the SEMF and 
QMC methods compared to the exact energy levels. The 
lighter areas indicate the QMC statistical errors. 
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FIG. 2. (Color online) SEMF, QMC (extrapolated esti- 
mate), and exact values of some observables for the yrast 
states of ^^Na of Fig. 1: Neutron (upper left) and proton (up- 
per right) occupations of the shells nlj normalized at 2j + 1, 
electric quadrupole moment (lower left), and magnetic dipole 
moment. The effective charges, orbital and spin g-factors 
have been chosen to be respectively e„ = 0.49e, gl^ = 0.036, 
g!^ — —3.875 for neutrons, and Cp = 1.29e, Qp = 1.03, 
gp = 5.586 for protons. 



Q2 and magnetic dipole /ii moments. The results for 
these observables as obtained with the SEMF and QMC 
methods are compiled in Fig. 2 for the nucleus ^^Na 



and compared to exact diagonalization. As above, the 
SEMF approach again provides a good approximation, 
the agreement with the exact values being roughly cor- 
rect for any angular momentum. Nevertheless, in con- 
trast to the case of the energy, no significant improvement 
of the SEMF results is observed after the imaginary-time 
propagation for these observables that do not commute 
with the Hamiltonian. This may point out the neces- 
sity to go beyond extrapolated estimates by using, for 
instance, the back-propagation techniques [19]. Finally, 
for all the states of Fig. 1, we have found exact expecta- 
tion values of the squared isospin. 

In summary, we have presented a new QMC formal- 
ism for the nuclear shell model. The method relies on a 
mean-field wave- function entangled by symmetry restora- 
tion to initiate, guide, and also constrain the Brownian 
motion in order to control sign/phase problems. The re- 
sults reported in this Letter demonstrate that the phase- 
less QMC approach can accurately reproduce the yrast 
spectroscopy of sd- and p/-shell nuclei with realistic effec- 
tive interactions. Electromagnetic transitions and /3 de- 
cay can also be considered via a specific mixed-estimator 
[20]. Finally, calculations on medium- mass neutron-rich 
and proton-rich nuclei are under development to investi- 
gate low- lying structure in the full pf — OgQ/2 shell that 
can only be addressed by conventional diagonalization 
methods through severe truncations of allowed configu- 
rations. We are also extending the phaseless QMC ap- 
proach for computing excited states of a given angular 
momentum in order to achieve a complete spectroscopy 
of nuclei. 

We gratefully thank R. Fresard and P. van Isacker for 
a careful reading of the manuscript. 
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